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ABSTRACT 

Context. A chemical scenario was proposed for photon-dominated regions (PDRs) according to which UV photons from nearby stars 
lead to the evaporation of very small grains (VSGs) and the production of gas-phase polycyclic aromatic hydrocarbons (PAHs). 
Aims. Our goal is to achieve better insight into the composition and evolution of evaporating very small grains (eVSGs) and PAHs 
through analyzing the infrared (IR) aliphatic and aromatic emission bands. 

Methods. We combined spectro-imagery in the near- and mid-IR to study the spatial evolution of the emission bands in the prototypical 
PDR NGC 7023. We used near-IR spectra obtained with the IRC instrument onboard AKARI to trace the evolution of the 3.3 /am and 
3.4 /am bands, which are associated with aromatic and aliphatic C-H bonds on PAHs. The spectral fitting involved an additional broad 
feature centered at 3.45 /am that is often referred to as the plateau. Mid-IR observations obtained with the IRS instrument onboard the 
Spitzer Space Telescope were used to distinguish the signatures of eVSGs and neutral and cationic PAHs. We correlated the spatial 
evolution of all these bands with the intensity of the UV field given in units of the Habing field Go to explore how their carriers are 
processed. 

Results. The intensity of the 3.45//m plateau shows an excellent correlation with that of the 3.3yum aromatic band (correlation 
coefficient R = 0.95) and a relatively poor correlation with the aliphatic 3.4 pm band (R=0.77). This indicates that the 3.45 pm feature 
is dominated by the emission from aromatic bonds. We show that the ratio of the 3.4 pm and 3.3 pm band intensity (13,4/133) decreases 
by a factor of 4 at the PDR interface from the more UV-shielded layers (Go ~ 150, 73 . 4 /Aj = 0.13) to the more exposed layers 
(Go > 1 X 10"^, 13,4/13,3 = 0.03). The intensity of the 3.3 pm band relative to the total neutral PAH intensity shows an overall increase 
with Go, associated with an increase of both the hardness of the UV field and the H abundance. In contrast, the intensity of the 3.4 pm 
band relative to the total neutral PAH intensity decreases with Go, showing that their carriers are actively destroyed by UV irradiation 
and are not efficiently regenerated. The transition region between the aliphatic and aromatic material is found to correspond spatially 
with the transition zone between neutral PAHs and eVSGs. 

Conclusions. We conclude that the photo-processing of eVSGs leads to the production of PAHs with attached aliphatic sidegroups 
that are revealed by the 3.4 pm emission band. Our analysis provides evidence for the presence of very small grains of mixed aromatic 
and aliphatic composition in PDRs. 

Key words. ISM: photon-dominated regions - ISM: individual objects: NGC 7023 - ISM: molecules 


1. Introduction 

A significant fraction of interstellar carbon (up to 20%, Joblin 
et al. 1992; Tielens 2005) is tied up in the carriers of the aro¬ 
matic infrared bands (AIBs). The most intense AIBs are ob¬ 
served at 3.3, 6.2, 7.7, 8.6, 11.3, and 12.1 pm and are generally 
attributed to stochastically heated polycyclic aromatic hydrocar¬ 
bons (PAHs, Leger & Puget 1984; Allamandola et al. 1985). 
This set of emission bands shows systematic variations in the 
relative intensity and central wavelength of the bands (Peeters 
et al. 2002), reflecting variations in the local physical conditions 
and a chemical evolution of their carriers. From the analysis of 
spectro-imagery data on several Galactic photon-dominated re¬ 
gions (PDRs), the variations of the mid-IR (5.5-15 pm) bands 
were attributed to the emission of three different populations: 
PAH cations and neutrals (PAH^ and PAH^) and evaporating 
very small grains (eVSGs, Rapacioli et al. 2005; Berne et al. 
2007; Pilleri et al. 2012). A chemical scenario has been pro¬ 
posed according to which eVSGs evaporate under UV irradia¬ 
tion, producing free-flying PAHs. Candidates such as PAH clus¬ 


ters (Rapacioli et al. 2006; Montillaud & Joblin 2014) or Fe-PAH 
complexes (Simon & Joblin 2009) have been proposed as mod¬ 
els for these eVSGs. Jones (2012) recently proposed a model for 
the evolution of carbonaceous dust in the interstellar medium in 
which the smallest particles are aromatized, whereas the largest 
ones remain H-rich and predominantly aliphatic. Such a scenario 
would reconcile the observation of both the aromatic emission 
bands in UV-irradiated regions and the 3.4pm CH aliphatic ab¬ 
sorption feature as a main feature of interstellar dust in galax¬ 
ies (Dartois et al. 2004). To develop this evolutionary scenario 
further, we report here a study of the 3.4pm emission band, 
which is observed as a minor band at the side of the main AIB at 
3.3 pm. Characteristic of aliphatic CH bonds, this band has been 
attributed to methyl sidegroups attached to PAHs (Jourdain de 
Muizon et al. 1990) or to superhydrogenated PAHs (Bernstein 
et al. 1996). The observed 3.4/3.3 pm intensity ratio ( 13 , 4 / 13 , 3 ) 
is found to decrease for an increasing intensity of the UV field 
(Geballe et al. 1989; Joblin et al. 1996; Sloan et al. 1997; Mori 
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Fig. 1: In color scale, the Spitzer IRAC 8 pm image of the NGC 7023 NW PDR (Werner et al. 2004). The dark lines demarcate the 
region in which AKARI-IRC data were obtained, the white rectangle represents the Spitzer-IRS field of view. The star represents 
the position of HD 200775. The insets represent the AKARI and IRS beam sizes. AKARI spectra between 2.8 and 3.9 pm are shown 
to the right. They correspond to the positions P1-P4 that are marked as black dots on the map. 


et al. 2014), refiecting the photo-destruction of the more fragile 
bonds linked to the 3.4 pm band. 

In this paper we present 2D spectroscopy of the PDR associ¬ 
ated with the reflection nebula NGC 7023 NW obtained with the 
AKARI and Spitzer space telescopes. The combination of these 
datasets allows us to study the spatial variations of the 3.3 and 
3.4 pm emissions due to aromatic and aliphatic bonds with those 
of the PAH^ and eVSG emissions. In Sect. 2 we present the ob¬ 
servations and the data reduction. In Sect. 3 we describe the data 
analysis tools and present the observational results. We discuss 
these results in Sect. 4. Finally, the conclusions are presented in 
Sect. 5. 


2. Observations 

2.1. NGC 7023 NW 

NGC 7023 is a reflection nebula in the Cepheus constellation il¬ 
luminated by the B2-5Ve binary star HD 200775 (Alecian et al. 
2008). The distance of HD 200775 to the Sun measured by 
Hipparcos is 430^^^^ pc, but a recent analysis revised the dis¬ 
tance to 320 + 51pc (Benisty et al. 2013). The nebula has been 
shaped by the star formation process, which led to the forma¬ 
tion of a cavity (Fuente et al. 1993). The interaction of the UV 
photons with the walls of the parent molecular cloud has pro¬ 
duced several PDRs to the northwest, south, and east of the star. 
The brightest PDR is found about 45" NW from the binary star 
(hereafter the NW PDR, or NGC 7023 NW). Pilleri et al. (2012) 
estimated the intensity of the radiation field at the PDR front 
to be Go ~ 2600 in units of the Habing field (Habing 1968). 
Previous observations of different gas lines have shown that this 
region hosts structures at different gas densities: ~ 100 cm“^ 

in the cavity (Berne & Tielens 2012), 10^-10^ cm“^ in the fila¬ 
ments that are observed in the mm (Fuente et al. 1996) and near- 
IR (Lemaire et al. 1996; Martini et al. 1997), and ~ 10"^ cm“^ in 
the molecular cloud (Gerin et al. 1998). 


2.2. Near- and mid-IR observations with Spitzer and AKARI 

The infrared camera (IRC, Onaka et al. 2007) onboard the 
AKARI space telescope (Murakami et al. 2007) was used to 
obtain near-IR (2.5-5 pm) spectro-imagery of NGC 7023 NW. 
The observations were performed as part of the open-time pro¬ 
gram NESID of the AKARI post-helium phase between October 
2008 and January 2010. They were carried out using the grism 
spectroscopy mode with the narrow Nh slit, which provides a 
spectral resolution of about 0.02yum and a spatial resolution 
of ~ 3" (Onaka et al. 2007). The slit position was slightly 
shifted for the different runs to make dithering observations. In 
each pointed observation, eight to nine exposure cycles in spec¬ 
troscopy were obtained together with one exposure in imaging 
in the N3 (3.2 pm) band, which gives the accurate position infor¬ 
mation. The on-source integration time for each observation is 
typically 360s. The data were processed with the latest version 
of the spectroscopic toolkit (version 20110301) for the phase 3 
observation^ The NGC 7023 NW observations consisted of 20 
stripes of length ~ 50", each with a slightly different position 
and orientation. The stripes were then combined in a mosaic us¬ 
ing the software montage^. The field of view of the final mosaic 
results from the overlap of the slit coverages of all the stripes 
and is displayed in Fig. 1. The near-IR features vary smoothly 
between nearby positions but are significantly different with dis¬ 
tance to the star, as shown by the four representative spectra in 
Fig. 1 that correspond to offsets of PI [-29", 24"]; P2 [-36", 
35"], P3 [-39", 38"], P4 [-43", 41"] relative to the position of 
HD 200775. 

We compare this dataset with the Spitzer-IRS observations of 
NGC 7023, which provides spectral imagery in the wavelength 
range [5.5-14] yum and covers the whole filamentary region of 
the NW PDR with a spatial resolution of 3.6"(see Fig. 1). This 
dataset has been described in detail in previous studies, for ex¬ 
ample, Werner et al. (2004), Berne et al. (2007), and Pilleri et al. 
( 2012 ). 


^ http://www.ir.isas.jaxa.jp/ASTR0-F/0bservation/ 
DataReduction/IRC/ 

^ http://montage.ipac.caltech.edu 
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Fig. 2: Summary of the observations and spectral analysis for positions PI to P4 (top to bottom row). First column: the observed 
spectrum (black) and its fit (red) for the AKARI observations. The fit of the 3.29 and 3.40 /im features are displayed in blue, and the 
3.45 yum plateau is shown in green. Second column: observed Spitzer-IRS spectrum (crosses), its fit (black) and the decomposition 
obtained with PAHTAT (red: eVSG, green: PAH^, blue: PAH^). 


3. Spectral analysis and spatial distributions 

3.1. Near-infrared 

The first column of Fig. 2 shows the AKARI-IRC spectra ob¬ 
served at the positions P1-P4. The spectra are dominated by the 
two emission bands that peak at 3.29 and 3.40/im, in addition 
to an underlying plateau. To fit these spectra we used two nar¬ 
row Gaussians for the 3.29 and 3.40yum features and a broader 
Gaussian centered at 3.45yum to account for the plateau. We let 
the central wavelength of the Gaussians vary by 0.01 yum to ac¬ 
count for small shifts and treated the band widths and intensities 
as free parameters. We simultaneously fit the observed spectrum 
using a linear combination of these features and a linear con¬ 
tinuum. The resulting full widths at half maximum of the bands 
were found to be ~ 0.03 yum for the 3.29 and 3.40 yum bands 
and ~ 0.1 yum for the 3.45 jum plateau. As shown in Fig. 2, this 
yields very good fits of the spectra, even in cases with a strong 
plateau and faint superposed 3.40yum band. 

The 3.29yum band is commonly attributed to the C-H in¬ 
plane stretching mode in PAHs, whereas the 3.40yum band and 
the plateau are more difficult to attribute unequivocally. Steglich 
et al. (2013) concluded that it is difficult from a spectroscopic 
point of view to favor the attribution of the 3.4yum band ei¬ 


ther to methyl side groups attached to PAHs or to superhydro- 
genated PAHs. Still, Joblin et al. (1996) showed that the evo¬ 
lution of the 3.4/3.3yum band ratio with the UV field can be 
explained by methylated PAHs. The 3.45yum plateau is most 
likely the superposition of a number of faint features includ¬ 
ing other aliphatic modes, hot band emission of aromatic CH 
(Barker et al. 1987) and combination bands of C-C and C-H 
aromatic modes (Allamandola et al. 1989). In Fig. 3 we show 
the pixel-to-pixel correlation of the integrated intensity of the 
3.45 yum plateau (Ipia) with that of the 3.29yum (I^ ^) and 3.40yum 
( 73 . 4 ) bands. Ipia correlates better with I 3 3 (correlation coeffi¬ 
cient R = 0.95) than with 73.4 (R = 0.77). This suggests that 
the plateau is dominated by the emission from aromatics (hot 
bands and combination bands). Since our goal is to study the 
behavior of the aromatic vs aliphatic evolution, we exclude the 
plateau from the analysis in the following. 

3.2. Mid-infrared 

To extract the emission of the different AIB carriers, we used 
the results of the PAHTAT procedure (Pilleri et al. 2012). Figure 
2 shows the mid-IR spectra and the corresponding fits for the 
four positions P1-P4. PAHTAT also allows estimating the value 
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Fig. 3: Pixel-to-pixel correlations of the integrated intensity of the 3.3yum band (I 33 ), 3.4yum band ( 73 . 4 ), and the 3.45yum plateau 
(Ipia)- The correlation coefficient R and the slope m of the linear fit to the points are also given. Ipia correlates better with the 
intensity of the aromatic 3.3 yum band than with the aliphatic 3.4yum band. This indicates that the plateau is dominated by emission 
from aromatic bonds. 


of Ay along the line of sight assuming that the AIB carriers and 
other dust populations are well mixed. From these Ay values, 
we reconstructed the deredenned mid-IR and near-IR spectra by 
using the extinction curve from Weingartner & Draine (2001) for 
an Ry value of 5.6 (Witt et al. 2006). These unreddened spectra 
are used in the following analysis. The 3.3 and 3.4yum bands 
are sufficiently close in wavelength for their ratio is to be not 
significantly affected by the extinction correction. 

Table 1 shows the relevant results of the fit for the posi¬ 
tions P1-P4, that is, the integrated intensity of the near-IR bands, 
of the PAH^, PAH^, eVSG, and of the total AIBs, as well as 
the fraction of carbon locked in eVSGs (/evsc)- In Pilleri et al. 
( 2012 ), we have shown that feysG can be derived from every 
pixel and then used as a probe for Go in the region of the PDR 
with significant eVSG emission. For NGC 7023 NW, this region 
starts at a distance of ~ 45" from HD 200775, corresponding 
to Go = 2600. In the more exposed layers of the PDR, we ob¬ 
tain an estimate of the UV field intensity by assuming geometri¬ 
cal dilution of the stellar radiation field, a stellar temperature of 
Teff = 15 000 K (Finkenzeller 1985), and a local extinction of 
Ay = 1.5 (Pilleri et al. 2012). The total map of Go we obtain is 
shown in Fig. 4. 


3.3. Spatial distributions 

Applied to the full IRS data cube, PAHTAT allows deriving the 
spatial distribution of the intensity of PAH^, PAH^, and eVSGs 
(top panel in Fig. 5). These results have been described in detail 
in Pilleri et al. (2012) and can be summarized as follows: PAH^ 
dominate the mid-IR emission between the illuminating star and 
the PDR front; PAH^ are the most abundant species at the PDR 
front as traced by the emission of the H 2 rotational lines; eVSGs 
dominate the emission slightly deeper inside the cloud; the tran¬ 
sition between PAH^ and eVSGs corresponds to the peak of the 
H 2 S(3) line, which is found at a magnitude of visual extinction, 
Av ~ 1 (Pilleri et al. 2012). 

Figures 5a,b show the integrated intensity maps of the 3.3 jum 
and 3.4yum bands. Both maps peak toward the H 2 S(3) filament 
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Fig. 4: Map of Go in NGC 7023 NW (color scale). The H 2 S(3) 
integrated intensity is shown for reference (linear steps of 1.7 x 
10“^ Wm“^sr“^). The values of Go between the star and the 
PDR were obtained assuming geometrical dilution of the stellar 
radiation field and a local extinction value of Ay = 1.5. Inside 
the PDR, we used the fraction of carbon locked in eVSG as a 
probe for Go (Pilleri et al. 2012). 


that delineates the PDR front. The 3.3 jam band is also detected 
in the cavity between the star and the PDR, whereas the 3.4yum 
band is detected only toward the PDR. Neither the 3.3yum nor 
the 3.4yum bands is detected in the more shielded layers of the 
PDR, farther away from HD 200775 than the H 2 filaments. 

In Fig. 5c,d we compare these maps with the spatial distri¬ 
bution of PAH^ and eVSGs intensities obtained with PAHTAT. 
A very good spatial correlation is found between I 3 3 and Ipah^ 
(Fig. 5c). In contrast, the 3.4yum band is more localized on the 
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Fig. 6: First column: Correlation between the observed I 33 IIPAHO and h.AlhAm as a function of Gq. Second column: Same plots for 
model maps of lip Am and hAl Ip Am computed using Eqs. 6 and 8 and the obtained values for IpAmllpAm (see text for details). 
Third column: Semi-empirical model that in addition includes the dependency with Go on the intensity at 3.3 and 3.4 yum relative 
to all PAH^ bands (see text). In all panels, the linear fit to the data points (red line) is shown as well as its slope m. 


Table 1: Results of the fit of the mid- and near-IR observations. 



Unit 

PI 

P2 

P3 

P4 

leVSG 

10 “^ Wm-2sr-‘ 

0.00 

4.81 

6.00 

4.01 

IpAlfi 

10 “^ Wm-2sr-‘ 

4.04 

11.4 

7.60 

2.63 

IpAH+ 

10 “^ Wm-2sr-‘ 

8.87 

2.13 

1.10 

0.73 

Iaib 

10 “^ Wm-2sr-‘ 

12.9 

18.3 

14.7 

7.36 

Ay 

mag 

<2 

15 

19 

23 

feVSG 

- 

0.00 

0.42 

0.58 

0.71 

Go 

- 

7000 

2600 

200 

150 

h .3 

10 -’Wm-2sr-‘ 

11.0 

15.8 

22.3 

6.47 

h.A 

lO-’Wm-^sr' 

0.31 

1.42 

2.34 

0.80 

^pki 

lO-’Wm-^sr' 

5.36 

13.5 

15.6 

4.80 

134/133 

- 

0.028 

0.090 

0.11 

0.13 


dense PDR with a good spatial correlation between I^a and levsG 
(Fig. 5d). Still, the 3.3yum band is observed inside the cavity. 


which indicates that there is a contribution from PAH^ to its 
emission, and/or that a fraction of PAH^ is present in the cavity 
(Fig.5e). The intensity at 3.3 yum decreases by a factor > 5 (from 
3x10“^ to 6x10“^ W m“^ sr“^) in regions where IpAmllpAm ^ 3 
(Fig. 5f). Since we are interested in the evolution of the 3.4yum 
band relative to the 3.3yum band, we restrain our analysis below 
to the PDR front and mask the pixels that have Ip Am 11PAm > 1 
and those where the 3.4 jum band is not detected. 

4. Discussion 

4.1. Variation of I 33 with Go 

The dependence of the 3.3 jum band intensity on Go results from 
three effects: change in the abundance of its carriers, PAH ion¬ 
ization since the 3.3 jum band is expected to be smaller for PAH^ 
relative to PAH^ (see a review on quantum-chemical calculations 
by Pauzat 2011), and excitation effects. We discuss these effects 
here by considering the I 3 3 jlpAm ratio. As shown in Fig. 6a, this 
ratio presents some dispersion, but with an overall increase with 
Go. 
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Fig. 5: Top: Spatial distribution of the 3.3yum (a) and 3.4yum 
(b) integrated intensity as observed by the IRC instrument 
(color scale). The contours represent the integrated intensity 
of the H 2 S(3) line at 9.7 yum (linear steps of 1.7 x 10“^ 
W m“^ sr“^), and the star indicates the position of HD 200775. 
Middle: Comparison of the integrated intensity of the bands 
at 3.3 and 3.4yum (color scale) with the spatial distribution of 
the PAH^ and eVSG intensity, respectively (contours steps of 
1.7 X 10“^ W m^ sr“^ starting at 1.7 x 10“^ W m^ sr“^). Bottom: 
Comparison of the integrated intensity of the 3.3yum band with 
that of PAH^ and with the ratio Ipah+ IIpah^ • In (f) the contours 
are shown in steps of 0.3, with increasing values while approach¬ 
ing the star from the PDR. 


where 63.3 is the average emissivity at 3.3yum per aromatic CH 
bond in PAHs along the line of sight, which depends on the 
molecular properties of their carriers and the excitation condi¬ 
tions, and A 3.3 is the column density of aromatic CH bonds. 

The intensity of the 3.3 yum band is the sum of the contribu¬ 
tion of the neutral and cationic PAHs. Thus, the ratio hjlIpAH^ 
can be written as 


^ 3.3 


IpAH^ ipAW 


h.i 


1 -h 


^3.3 

^3.3 


3.3 


'^3.3 


^PAW €p 


PAH 


1 -r 


N: 


3.3 ^3.3 
3.3 3.3 


( 2 ) 


where 3 ^, 633 ^ are the emitted intensity and average emissiv¬ 
ity of the 3.3 yum band due to PAH^ and PAH^, and A 3 3 ^ is the 
column density of aromatic CH bonds of neutral and ionized 
species. We assume that the ionization fraction is the same for 
the carriers of the 3.3yum band and for the whole PAH popula¬ 
tion, that is, A 33 /A 33 = NpAH+INpAH^ (Montillaud et al. 2013, 
see discussion below). From this we obtain 


h.3 j I. NpAH+ 
-— A 3 311 -f- 

hAW \ ^PAH^ 

where we have defined 



( 3 ) 


^ 3.3 


*^3.3 


A^ 


4aH 


n - 

P33 - 

*=3.3 


( 4 ) 

( 5 ) 


^ 3.3 represents the fraction of the bolometric PAH emission 
in the 3.3yum band for the PAH^ population and is therefore an 
indicator of the PAH size distribution. Pech et al. (2002) have 
modeled the IR emission spectrum of a distribution of PAHs 
under irradiation by UV photons and derived values for ^ 3.3 
of about 0.01 - 0.02 depending on the assumed size distribu¬ 
tion. The 773.3 parameter represents the emissivity ratio for PAH 
cations and neutrals. It is expected to depend on the excitation 
conditions since PAH cations (open-shell species) have more 
electronic bands at low energy. However, in a relatively hard UV 
field, the excitation will be dominated by the strong UV n ^ n* 
and (T ^ (T* type bands, which are present in all charge states 
and are almost equal across them (Mulas et al. 2006). The 773 3 
parameter will then depend mainly on the relative IR absorp¬ 
tion cross-section at 3.3yum of PAH cations relative to neutrals. 
Values in the range 0.2-0.5 have been calculated for PAH sizes 
between 20 and 60 C atoms (Malloci et al. 2007). 

The total emissivity for a neutral and a cation with similar 
size is approximately the same because of conservation of en¬ 
ergy and the fact that their UV absorption cross-sections are in 
first approximation similar, as explained above. Therefore, we 
can use the ratio of the observed intensity Ipah+ I^paw as a probe 
for the column density ratio of ionized to neutral PAHs, and 
rewrite Eq. 3 using only the observed intensities and the molec¬ 
ular parameters ^ 3.3 and 773 . 3 : 


The observed intensity of the 3.3 yum band can be written as 
a function of the column density A 3.3 of its carriers and Gq: 

The term in parenthesis is related to the contribution of 
A .3 = 670 A 3 . 363 . 3 , (1) the ionization fraction of PAHs along the line of sight to the 


h3 

PAH^ 


= k33 


4 IpAH+ 

1 + - - P33 

^PAH^ 


( 6 ) 
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observed 133 /IpAH^ ratio. Using the maps we obtained with 
PAHTAT, we derived the values of Go and Ipah+ /Ipah^ for each 
pixel (Figs. 4 and 5f, respectively). Then, we used these values 
to compute a model map of h.^jlpAH^ using Eq. 6 with con¬ 
stant values of ^ 3.3 and 773 . 3 . This semi-empirical model (here¬ 
after model 1, Fig. 6 b) cannot reproduce the observed trend in 
h.^lIpAH^- As an example, Fig. 6 b shows the results for model 1 
using ^ 3.3 = 0.2 and 773.3 = 0.5. These values give the best fit to 
the data within the physically meaningful range for ^ 3.3 and 773 . 3 . 
This shows that the observed variations of h.^jlpAH^ cannot be 
accounted for by a variation in the mean charge of PAHs alone. 

To obtain a better fit of the observed increase of I 33 IIpah^ 
with Go, we assumed that ^ 3.3 depends on Go with a linear rela¬ 
tionship: ^ 3.3 = 0.02 -r /:* 3 X Go- We fit the slope of this correla¬ 
tion by varying 3 . The best fit is obtained using ^3 3 = 8 x 10“^ 
(hereafter model 2, Fig. 6 c). The dependency of ^ 3.3 on Go could 
be due to an increase of the abundance of the 3.3 7 /m carrier rel¬ 
ative to all neutral PAHs, or to an increase of the emissivity at 
3.3 7 /m relative to the total emissivity of PAH^. As discussed 
below, the first effect is likely to involve a chemical evolution of 
the PAH population, whereas the second is related to excitation 
conditions. 

Analyzing infrared observations obtained by Spitzer and 
Herschel, Berne & Tielens (2012) showed that there is no global 
variation of the PAH abundance in the region covered by our ob¬ 
servations. Montillaud et al. (2013) modeled the chemical evolu¬ 
tion of PAHs in NGC 7023 NW. They showed that PAHs of small 
sizes (with a number of carbon atoms Nc < 50) are fully de¬ 
hydrogenated, whereas larger PAHs are normally hydrogenated 
and possibly superhydrogenated for the largest sizes (Nc > 90). 
The ionization fraction is found to be rather similar from one 
size to the other. It is interesting to note that across the PDR, 
intermediate-sized PAHs (Nc ~ 60-70) can experience strong 
variations in their hydrogenation state because of the competi¬ 
tion between photodissociation by UV photons and recombina¬ 
tion with H atoms. For instance, C 66 H 20 is predicted to be fully 
hydrogenated at the border of the PDR, but fully dehydrogenated 
deeper inside the cloud at a distance from 5 to 10" from the 
PDR front. Here, the abundance of H relative to H 2 decreases 
faster than the attenuation of UV photons, resulting in fully de¬ 
hydrogenated species. This effect could explain the increase in 
the 3.3 7 /m band intensity with increasing Gq. Another possibil¬ 
ity is that the emissivity at 3.3 7 /m increases with Go as a result 
of the higher average internal temperatures of PAHs with higher 
UV field. Two effects may be invoked: in the most exposed lay¬ 
ers of the PDR, the UV radiation field is becoming harder be¬ 
cause of decreased extinction by dust (Rapacioli et al. 2006). 
Secondly, when the UV flux increases, the probability for mul¬ 
tiple photon events increases. A PAH that has absorbed a UV 
photon can only relax part of its internal energy before absorb¬ 
ing another photon, and it is therefore heated to higher temper¬ 
atures. However, whereas multiple photon events are important 
for studying photo-dissociation (Montillaud et al. 2013), they re¬ 
main rare events in NGC 7023 NW and do not contribute signif¬ 
icantly to the mid-IR emission. In addition, Witt et al. (2006) 
showed that in the first layers of the PDR there is a significant 
contribution of high-energy UV-photons (hv > 10.3 eV) that are 
suppressed in the more protected layers (Ay > 2 ). The variation 
of the hardness of the UV field across the PDR provides another 
explanation for the increase of the intensity of the 3.3 7 /m band 
with Go. Quantifying the importance of this process compared 
to the hydrogenation effect described above would require a full 
chemical and photo-physical model, which is not the purpose of 
this paper. We conclude that the increase of h.^lIpAH^ with Go is 


due to a change in the hardness of the UV field combined with 
an increased abundance of aromatic CH bonds at the border of 
the PDR, both effects being relatively moderate. 

4.2. Photo-destruction of aliphatics 

Similar to the 3.3 7 /m band, the intensity of the 3.4 7 /m band rel¬ 
ative to the total intensity of the PAH^ bands (hAlIpAH^) also 
presents some dispersion. However, it shows an overall decrease 
with Go (Fig. 6 d). In analogy with the 3.3 7 /m case, the observed 
intensity of the 3.4 7 /m band can be written as 

I 3 A = G 0 N 3 A 63 A, (7) 

where 63.4 and V3.4 are the average emissivity and column den¬ 
sity per aliphatic CH bond along the line of sight. Equation 6 can 
be rewritten for the 3.4 7 /m band as 


73.4 

hAH^ 


= k3A X 



hAH+ 

IpAH^ 


XP3A 


where we have defined 


( 8 ) 


^ 3.4 = 


N^ 

3.4 3.4 


/O 

^3.4 


4ah 4ah 


n - 3.4 

P3A - 

*=3.4 


(9) 

( 10 ) 


To obtain an estimate of the typical values for 773 4, we used 
the NASA-Ames PAH IR Spectroscopic database^ (hereafter 
the NASA-Ames PAH database, Boersma et al. 2014) and ob¬ 
tained 773.4 > 1, depending on the size of the PAH being consid¬ 
ered. Since 773.4 > 773.3, the term in parenthesis is larger for the 
3.4 7 /m than for the 3.3 7 /m band, implying that the charge state of 
PAHs modifies the ratio hAlhAH^ more significantly than does 
h.3lIpAH^- 

Assuming 773.4 = 1, the best fit to the data using our model 1 
is obtained using ^ 3.4 = 0 . 002 , which allows us to reproduce 
the average values of I 3 a/Ipah^ (Fig. 6 e). ^ 3.4 is found to be 
a factor of 10 lower than ^ 3.3 because of the lower abundance 
of aliphatic CH groups compared to the aromatic CH groups. 
The relative abundance of aliphatic vs aromatic groups we derive 
agrees with previous studies, including the recent work of Li & 
Draine (2012). The observed negative slope of the I 3 aIIpah^ vs 
G o correlation (Fig. 6 d) is not reproduced by our model 1 (Eq. 8 
and Fig. 6 e). As for the 3.3 7 /m band, we assume a dependency 
of ^ 3.4 with Gq: ^ 3.4 = 0.002 -\-k*^x Gq. The observations can be 
best reproduced by ^34 = -1.8 x 10“^ (Fig. 6 e). This negative 
term mainly reflects the destruction of the carriers of the 3.4pm 
band with increasing Gq. This is an effective value that reflects 
both the increase in intensity and hardness of the UV field while 
reaching the edge of the PDR. 


4.3. Variation of I3 A I h. 3 with Go 

The ratio / 3 . 47 / 3.3 can be used as a tracer of the aliphatic to 
aromatic content in PAHs if we assume that both bands arise 
from the same carriers. This ratio is low in the cavity (~ 0.03 
at PI) and increases monotonically toward the PDR, reaching 
its highest value (~ 0.15) slightly behind the H 2 S(3) and the 

^ http://WWW.astrochem.org/pahdb/ 
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Fig. 7: Map of the 3AI33 yum ratio (color scale) with the contours of the H 2 S(3) line {a), and of the ratio levscIhAH = 
leVSclilpAH^ + hAH+) (^)- The density plot of the l 3 Alh .3 band ratios vs Go is also shown (c). In red, we report the linear re¬ 
gression of the data points. Only the pixels for which I 3 3 > 1.0 x 10“^W ni“^ sr“^ have been used to build in the analysis of the 
3A/3.3 yum ratio. 


PAH^ emission (Fig. 7a). The highest gradient is found at the 
PDR front where in only ~ 10" the ratio increases by a fac¬ 
tor of ~ 3. This shows that photo-chemical processing of the 
carriers of the 3.4yum band is taking place in this region. We 
note that 10 " corresponds to the angular width of the filament 
traced by the H 2 rotational emission. Interestingly, this is also 
the evaporation zone of eVSGs into gas-phase PAHs. We have 
used the levs d Ip ah - hvsG/ilpAH^ + lpAH+) ratio as an indi¬ 
cator of the destruction of eVSGs into PAHs by UV photons. 
The relative spatial distributions suggest that the destruction of 
eVSGs is then followed by that of the 3.4yum carriers. This sug¬ 
gests that the evaporation of eVSGs leads to the production of 
PAHs with aliphatic groups. 

Figure 7c shows the pixel-to-pixel correlation between 
hAlh .3 and log(Go). Dividing each side of Eq. ( 8 ) by Eq. ( 6 ), 
we obtain 


where Dion is the ratio of the terms in parenthesis in each equa¬ 
tion and represents the variation of l 3 A/h .3 hue to ionization. 
Because the intensity of PAH^ relative to that of PAH^ is < 0.4, 
Dion is very close to 1. More specifically, Dion is about 10% 
higher in the more exposed regions in which Ipah+ /Ipah^ = 0-4 
compared to the region with no PAH^ emission, thus it cannot 
explain the trend shown in Fig. 7c. For PAH sizes larger than 
~ 50 carbon atoms, 63 4/63 3 can be also considered constant (cf. 
NAS A-Ames PAH database), and therefore the observed trend 
can be mainly attributed to a decrease in the abundance of the 
aliphatic CH bonds relative to aromatic CH bonds. 

The trend we observe is similar to that reported in Joblin 
et al. (1996) and analyzed by the authors using a simple photo¬ 
chemical model that describes the evolution of methylated 
PAHs. The authors showed that the decrease of the 3.4yum band 
occurs because the reactions that could reconstruct the aliphatic 
side-groups in the gas-phase are highly improbable. 

5. Conclusions 

Thanks to the high spatial resolution of observations made in 
the near- and mid-IR, we probed the evolution of aliphatic and 


aromatic CH groups in the spatially resolved PDR NGC 7023 
NW. We showed that the intensity of the 3.3yum band relative 
to the total PAH emission increases with Go, while the rela¬ 
tive contribution of the 3.4yum band decreases with Gq. In the 
most exposed layers of the PDR, the UV radiation field is harder 
and high-energy UV photons can excite PAHs to higher tem¬ 
peratures, leading thus to an increase of the fraction of energy 
emitted at 3.3yum. The higher abundance of H atoms in this re¬ 
gion can also increase the abundance of CH aromatic bonds. On 
the other hand, the higher fiux of UV photons leads to efficient 
destruction of the more fragile aliphatic CH bonds attached to 
PAHs, which cannot be efficiently reformed in these regions. In 
addition, we showed that the observed 3.45yum plateau is domi¬ 
nated by the emission from aromatic bonds. 

We find that the change in the aliphatic vs aromatic com¬ 
position is particularly important along the filament that delin¬ 
eates the PDR. In particular, the l 3 A/h .3 ratio peaks close to 
the region where eVSGs are photo-evaporated into PAHs. This 
suggests that the processing of eVSGs leads to PAHs with at¬ 
tached aliphatic sidegroups, which provides further insights into 
the formation and evolution of these species in astrophysical en¬ 
vironments. (Very) small grains of mixed aromatic and aliphatic 
composition are known to be present in circumstellar environ¬ 
ments (Kwok et al. 2001; Goto et al. 2003; Sloan et al. 2007), 
and they have been considered to be a major component in the 
dust model by Jones (2012). We provided here the first obser¬ 
vational insights for the presence of such grains in interstellar 
clouds. 
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